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Abstract. The heavy quark system is an excellent probe to learn about the 
' QCD dynamics at finite density. First, we discuss the properties of the J/ip and 

^ ' D meson at finite nucleon density. We discuss why their properties should change 

' at finite density and then introduce an exact QCD relation among these hadron 

properties and the energy momentum tensor of the medium. Second, we discuss 
attempts to calculate charmonium-hadron total cross section using effective 
(N , hadronic models and perturbative QCD. We emphasize a recent calculation, where 

the cross section is derived using QCD factorization theorem. We conclude by 
' discussing some challenges for SIS 200. 

> 

(N 

. 1. J/^p and D mesons at finite density 

o . 

_ In the heavy quark system, the heavy quark mass provides a natural normahzation 

' scale 11, which makes the perturbative QCD approach possible even for calculating 

some bound state properties of the heavy quarks. This implies that the heavy quark 
system is an excellent probe to learn about the QCD dynamics at finite density and/or 
temperature. Here, we will highlight some calculations which explicitly demonstrate 
^ ' the relations between the changes of the properties of the heavy quark system and the 

pH I QCD dynamics at finite density. 



X 



1.1. J/ijj mass at finite density 



5_j ■ Let us first start with the J/t/j meson mass at finite density. The interaction between 

I the J/tjj and a hadron h can be perturbatively generated by multiple gluon exchanges. 

This implies that the mass shift of the J/ijj in nuclear medium can be estimated if the 
gluon distribution is known in nuclear medium Q. In fact, in the heavy quark mass 
limit TOc — > oo, the mass shift can be calculated exactly to leading order in QCD |^ 
and is given as 

Ato,//^ = Ci X A (S^) matter, (1) 

where ci is a calculable constant and the matrix element is taken with respect to the 
nuclear matter. To the leading order in density, the matrix element is known from 
the trace anomaly relation and the gluon distribution function of the nucleon |^ . For 
a more realistic charm quark mass, there are some model dependence on how to treat 
the bound state part. So far, there are potential model approaches ^ and QCD 
sum rule approaches ^. In all cases, similar relation to Eq. (|l|) holds with slightly 
different expressions for ci, such that the mass shift at nuclear matter density ranges 
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from —4 to —7 MeV. The mass shift is small at nuclear matter density and one doubts 
whether such a small mass shift can be observed. However, one can turn the argument 
around and claim that if any mass shift for J/ip is observed, it will tell us about the 
changes of gluon field configuration at finite density through Eq. (|l|). 

1.2. D meson mass at finite density 

In the D meson, the heavy quark acts as a source for the light quark, which surrounds 
the heavy quark and probes the QCD vacuum. This is the basic picture of the 
constituent quark in the heavy-light meson system within the formulation based on 
the heavy quark symmetry Q. Therefore, if the vacuum properties are changed at 
finite temperature or density, the light quark should be sensitively affected and be 
reflected in the changes of the D meson properties in matter. For the D meson, there 
does not exist a formal limit where one can calculate the mass change in matter. 
Nevertheless, model calculations suggest ^ ^ that it is dominantly related to the 
change in the chiral condensate. In fact, in both the quark-meson coupling model |^ 
and in the QCD sum rule approach ||l^, the average mass shift of the mesons at 
normal nuclear matter density were found to be around —50 McV. 

The D meson mass shift has several importance in the J/^ suppression 
phenomena in relativistic heavy ion collisions (RHIC). First, it leads to subthreshold 
production and the change of the DD threshold for the J/^ system in matter JTT| . 
Second, a change of the averaged mass might lead to level crossings of the DD 
threshold with the ■;/;' and x mass |l^ , which leads to a step- wise J/ip suppression 
Q m. Unfortunately, for the D meson mass shift, there are large model dependence 
and nontrivial splitting between the Z?+ and D^ mesons p^ . 

1.3. Exact sum rule from QCD 

The D meson mass and the DD threshold have important phenomenological 
consequences. Unfortunately, only model-dependent calculation could be made. Here, 
we will introduce an exact QCD equality, which can be used to link the DD threshold 
to J /ij) suppression. The exact sum rules in QCD at finite temperature and/or density 
were first introduced to the correlation functions between light hadrons and the 
derivation goes as follows. Let us consider the correlation function between two meson 
currents and its dispersion relation. 



AH(g2,r,p) = I J d^xe"'^A{qrq{x)qrq{0)) 



1 f , AlmH(s) 

Here, we have assumed = (iQ, 0,0,0) and A(-) = (•)mcdium - (•)vacuum means 
the difference between the medium and the vacuum expectation values. F is a Dirac 
gamma matrix that will be determined by what meson we want to study. At ^ oo, 
the real part of Eq. (||) can be determined from the operator product expansion (OPE) 
in QCD, 



Here, 7„ (7,7) is the anomalous dimension of the operator (current). Equation (|^) is 
an asymptotic expansion in l/log(Q-^) and l/Q^. For 7„ = 7j = 0, the leading term 
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J/\|/ v' DD 
Figure 1. Charmonium spectral density. 



is proportional to and can be related to the first moment of the spectral density 



as 



ds-AImn(s) 

TT 



C4A(04), 



(4) 



where O4 is a dimension 4 operator and C4 its corresponding Wilson coefficient. 

One can generalize this formula for the heavy (charm) vector current 07^0 at rest 
and obtain the following constraint: 



/ ds-AImn(.) = - J-A(^G2) + ^,^A(roo) 
J TT 12 n 16 + 3nf 



(5) 



where Tf^^, is the traceless energy momentum tensor and (^G^) is the gluon 
condensate, which is related to the trace of the energy momentum tensor. Lattice 
gauge theory provides the temperature dependence of the operators both below and 
above the phase transition temperature. From this, one can investigate what kind 
of changes in the spectral density, which in the vacuum has a generic form given in 
Fig. ^, is consistent with the constraint equation. 

Now let us look at the changes at finite density. Using the linear density 
approximation A(-) = 2m (p\ ' '^^'^ 

1 . ,Q! 



6.25 X 10-^ X X GeV* 



(6) 



16 + 3n/ 



A(Too} « 3 X 10"'' X X GeV*, 



where x = p/pn.m. with pn.m. — 0.17/(fm)^. In comparison with this, when neglecting 
the contribution from the ip' , the spectral density in Fig. |l| can be parameterized as 



TT 



where = 



Imn(s) 

r(J/V' 



e+e-) 



1 



(1 



ds-AImn(s) 

TT 



Att^ ' Att' 
0.16 GeV'^. Hence, we have 
1 



J/i' 



d,TT' 



(1 + ^)A(^^). 



(7) 



(8) 



Now, since the mass shift ts.mj is very small, the change in the first term is coming 
from the change in /y , which is related to the change of the coupling of the J j t/i to the 
electromagnetic current. This is precisely related to the amount of J/V' suppression 
in matter. At linear density approximation, as can be seen from Eq. (6), the expected 
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change in the OPE is of 10^"' GeV' in the nuclear matter and can be neglected 
compared to a fractional change in the phenomenological side of Eq. (^). Therefore, 
in this case, the suppression in fy is related to the decrease of Sq. In an order of 
magnitude estimate, we arrive at 

where Afy/fy quantifies the amount of J/^p suppression seen through the dilepton 

signal. Assuming that Sq/' = mo + rnjj, Eq. ^ implies that 10% suppression is 
related to 28 MeV mass decrease of the D meson in nuclear matter. 

The relation in Eq. (||) can be used in the other ways. Namely, from the 
measurement of the changes in the mass of the D mesons and the suppression of 
the J /ip peak, one is able to learn about the changes in the energy momentum tensor 
and gluon condensate at finite temperature or density, where the measurement would 
be made. 



2. J/V'-hadron cross section 

J/il) suppression ||T^ seems to be one of the most promising signal for QGP formation 
in RHIC. Indeed the recent data at CERN show an anomalous suppression of J 
formation, which seems to be a consequence of QGP formation iQ. However, before 
coming to a concrete conclusion, one has to estimate the amount of J /Tp suppression 
due to other non-QGP mechanism, one of which is the hadronic final state interactions. 
Indeed, there are model calculations that show large suppression by hadronic final 
state interactions |2^, |2^ . The typical values for J/V' + hadron dissociation cross 
sections one uses for these calculations are (Jj/xp+N ~ 5 mb and (Jj/^+comover ^ 1 nib. 
However, there is no direct experimental data from which we can determine or confirm 
these cross sections. In fact, the existing model calculations vary greatly in their 
energy dependence and magnitude near the threshold. In Fig. |[ we show the latest 
typical results of model calculations based on perturbative QCD [p2| , |2^ , ps] , |26|| , 
the meson exchange model |^ and the quark exchange model [^8[ . To understand 
the discrepancies and to obtain a consistent result, it is crucial to probe each model 
calculations further so as to spell out the limitations and range of validity of each 
model calculations. Here, we will have a closer look at the perturbative QCD result. 



2.1. Perturbative QCD methods 

More than two decades ago, Peskin and Bhanot and Peskin |2^ showed that 
one could apply perturbative Quantum Chromodynamics (pQCD) to calculate the 
interactions between the bound state of heavy quarks and the light hadrons. Such 
calculation was feasible because in the large quark mass limit one could consistently 
obtain the leading order OPE of the correlation function between two heavy meson 
currents in the light hadron state. The justification of such calculation stems from the 
fact that there exist two relevant scales in the bound state |^ in the large quark mass 
limit, namely the binding energy, which becomes Coulomb- like and scales as mg'^ ^ and 
the typical momentum scale of the bound state, which scales like mg^, where m is 
the heavy quark mass and g the quark-gluon coupling constant. Hence, taking the 
separation scale of the OPE to be the binding energy, it is consistent to take into 
account the bound state property, which is generated by the typical momentum scale 
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Figure 2. u | ^ cross section from meson excliange model |27| , quark 
exchange model p8[ | and perturbative QCD [ p3| calculations. 



of the bound state, into the process-dependent Wilson coefhcient. The result obtained 
by Peskin and Bhanot was derived within the OPE [^2[ . Recently, we have derived 
anew the leading order pQCD result using the QCD factorization theorem The 
equivalence between the OPE and the factorization approach is well established in 
the deep inelastic scattering and Drell-Yan processes [|o|. Our work shows that the 
equivalence between the two approaches are also true for bound state scattering to 
leading order in QCD. Similar approach has been used by one of us (SHL) to estimate 
the dissociation cross section of the J /il) at finite temperature The factorization 
formula also provides a manageable starting point to calculate higher twist gluonic 
effects ^, which should be nontrivial for the J/-0-hadron scattering. Here, we will 
sketch the derivation . 

We refer to a bound state of heavy quark and its own antiquark as $. According 
to the factorization formula, the total scattering cross section of $ with a hadron h 
can be written as 

c<t,h{v) = / dxa^g{xiy)g{x), (10) 
Jo 

with v = p ■ q/M^, where p (q) is the momentum of the hadron (<i>) and is the $ 
mass. Here, a^g is the perturbative $-gluon scattering cross section and g{x) is the 
leading twist gluon distribution function within the hadron. The separation scale is 
taken to be the binding energy of the bound state. Hence, the bound state properties 
have to be taken into account in cr$g. This can be accomplished by introducing the 
Bethe-Salpeter (BS) amplitude r(pi, — P2), which satisfies [|3| 

X r^(Pi + k, -P2 + k)A{-p2 + k)"f^, 

(11) 

where Ccoior — {N^ ~ l)/(2^c) with the number of color Nc- The trivial color indices 
have been suppressed and pi (— P2) is the four-momentum of the heavy quark (anti- 
quark). 



T^{pi,-p2)=tC,olor / — -45V(fc)7''A(pi fc) 
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Figure 3. Scattering processes of "I> and gluon. 



We write q — pi + P2 and p — (pi — P2)/2 and work in the $ rest frame and 
introduce the binding energy e, such that qo — Af$ = 2m + £ (e < 0). Then, in the 
non-relativistic limit, the BS amphtude reduces to the following form: 



r a , 1 , ^ / If ^l±7o . 1 - 7o 



(12) 

and the corresponding BS equation becomes the non-relativistic Schrodinger equation 
for the Coulombic bound state, so that '0(p) is the normalized wave function for the 
bound state. 

With the BS amplitude defined as Eq. we now obtain the <i>-gluon scattering 
amplitude with the processes depicted in Fig. ^ The scattering matrix elements can 
be obtained by substituting the BS amplitude for the vertex. To obtain the leading 
order result, we should pay attention to the counting scheme. First the binding 
energy £q = m [iVc.g^/(167r)] is 0{mg^). Combined with the energy conservation 
q + k = pi + p2 in the non-relativistic limit, this implies |pi| ^ |p2| 0{mg^) and 
A;" = |k| 0{mg^). Within this counting scheme, the final answer reads 



X "(Pi)^— 7 -^—^ "(P2)- (13) 
The scattering cross section a^g can now be obtained from 

With the amplitude (^3|) we obtain 

Substituting Eq. ([l^) into Eq. (|l^) and using the known gluon distribution function 
in the pion, we obtain the dashed line in Fig. ^ 

2. 2. Corrections to leading order formula 

The counting of this formalism implies that for the charmonium system = 0.84 
and the separation scale Eq = 0.78 GeV. Hence, both the and higher twist effects 
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Figure 4. Scattering cross sections aj/^_g and fi/,'(2S)-g- The solid lines are 
obtained in the non-relativistic limit and the dashed lines are with the relativistic 
correction. 



should be important. However, calculating higher corrections in this formalism is 
intricately related to the relativistic corrections including 1/m corrections, which is a 
formidable future task. 

Here, we will only look at the simple but important relativistic correction coming 
from the relativistic calculation of the phase space integral in Eq. (|lj) . The difference 
after a full calculation is shown in Fig. ^ for cr^g, which shows larger suppression of the 
cross sections at higher energy. However, we find that this relativistic correction has 
almost no effect on the cross section cr$/i. This follows because CT$g is highly peaked 
at small energy region and the gluon distribution function increases towards small x. 
Through Eq. (^, this implies that (T$/j, at all energies, are dominated by the low 
energy behavior of cr$g, which is reliably calculated in the non-relativistic limit. This 
implies that the full relativistic corrections to the scattering amplitude cr$g would have 
little effects on a^h- 

The remaining important and interesting correction is the higher twist-effect, 
which contributes as A/e'^ in our formalism. Since A is related to some hadronic scale 
of the hadron h, it should be of the order itself, hence should be non-negligible 
. It is also important to extend this calculation to investigate the absortion cross 
section for the x states. This will influence the amount of J/ip production coming 
from the decay of the x's in a p-A or A-A reaction js^. 

3. Conclusion and challenges for SIS200 

As we have discussed, observation of the changes in the properties of the heavy quark 
system will provide us useful information about the QCD dynamics both in the vacuum 
and in the medium. Important experimental observation would be mass shifts of the 
charmonium and open-charm system in nuclear medium. Also important would be the 
quantitative observation of the suppression of the J/ip and ip' signals in the medium 
through the dileptons. 

p-A experiments to probe the J/'ip in nuclear matter would also provide useful 
informations on the J/'i/'-nucleon cross section at low energy. This is precisely where 
there is huge discrepancy between existing theoretical predictions. Experimental 
investigation will provide useful information on the higher twist effects in the pQCD 
result for the J/i/'-nucleon cross section. Here, the higher twist effect comes solely 
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from gluon operators and hence this in turn would be a unique place to learn about 
the higher twist gluon effect of the nucleon and nuclear matter. 
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